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Abstract 

In 1975, a modification was formulated to Besseling' s Subvolume Method to 
allow it to use multilinear stress-strain curves which are temperature 
dependent to perform cyclic thermoplasticity analyses. This method 
automotically reproduces certain aspects of real material behavior important 
in the analysis of Aircraft Gas Turbine Engine (AGTE) components. These 
include the Bauschinger effect, cross-hardening, and memory. This 
constitutive equation has been implemented in a finite element computer 
program called CYANIDE which has been in production usage since 1977. 
Subsequently, classical time dependent plasticity (creep) was added to the 
program. Since its inception, this program has been assessed against 
laboratory and component testing and engine experience. The ability of this 
program to simulate AGTE material response characteristics has been verified 
by this experience and its utility in providing data for life analyses has 
been demonstrated. In this area of life analysis, the multiaxial 
thermoplasticity capabilities of the method have proved a match for the actual 
AGTE life experience. This paper will explore the multiaxial, 
variable-temperature nature of the method and show examples demonstrating its 
utility. 
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BESSELING'S MATHEKATICAL MODEL 


The relation between the devlatorlc stresses and the devlatorlc 
strains Is given by 

- 2GCe^j - ejj) (1) 

where 

Is the devlatorlc stress tensor 
e^j Is the total devlatorlc strain tensor 
e^j Is the plastic strain tensor 
G is the shear icodulus 


The yield strain, P, is given by the plastic potential function 
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The Incremental plastic strains are given by 




provided that 
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The incremental stress-strain relation Is 
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The new vield strain, e^^ + is determined from 


6g - 2(e^j - “ «®Ij> " 0 


( 6 ) 


Bessellng then Introduced the concept of elastic-perfectly plastic 
subvolumes. The elastic potential, of the subvolume of density p 
after prior plastic flow is given by 
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where the e^^^^ are the plastic strains due to Ideal plastic yielding. 
If this subvolume constitutes the fraction i|/ of the volume element 
dV, Its contribution to the total elastic potential of dV Is 




(8) 


If k subvolumes of the volume dV have exceeded their critical value of 
elastic potential and undergone plastic flow, the total elastic 
potential Is given by 


p*dV 
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Now, the devlatorlc stress tensor Is given by 
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After yielding, the elasticity llolt of subvolume k Is given by 


- »ijk> - '■k - 0 


(11) 


The subvolume Incremental plastic strains are given by 
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provided that 
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The Incremental stress-strain relations are 
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DEVELOPMENT OF NON I SOTHERMAL CAPABILITY 


The equation relating the stresses and the subvolume strains, 
Equation (10), can be rewritten to give 
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Now these stresses must be the same as the stresses given by Equation (1) . 
Therefore, the two right-hand sides can be equated. When this 
Is done, we get 
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( 16 ) 
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which gives a relationship between the subvolune plastic strains and 
the total plastic strains. 


Squaring both sides of Equation (16) and multiplying by 2/3, we 


get 
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Now 


e 

P 




(18) 


Therefore 
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This gives a relationship between the total effective plastic strain 
and the subvolume effective plastic strains. 


The following ratio can be formed between a subvolume effective 
plastic strain and the total effective plastic strain: 


or 
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By taking the square root of both sides, we obtain 


'ijn 
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( 22 ) 


This gives a means of determining the subvolume plastic strains from 
the total plastic strains if the effective plastic strains are known. 


This then provides the tools to convert Bessel Ing's Isothermal 
theory Into a nonisothermal theory. We note that for variable 
temperature problems g and g^^ will be functions of both strain and 
temperature. 


9 *= g 

Sk ' Sk<'lJ-'^> 

These functions can be specified by defining temperature dependent 
stress-strain curves. 


( 23 ) 
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For Incremental loading Including temperature changes, the 
change in the plastic potential function Is given by 
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( 25 ) 


There are three possible conditions that can occur due to this 

load Increment and these are determined by the value of this differential. 


For loading beyond the present yield surface 
dg > 0 


* -If > 0 


( 26 ) 
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For the loading to place the point on the new yield 
surface 

dg = 0 (28) 
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For the point to unload back into the elastic range 
dg <0 

3g de.. + a|2_ dT< 0 

3e.. 3T 

•J 

These last two conditions are used to accomodate 
temperature variations. The solution to any load condition, 
(N-1), is arrived at when 
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T = Constant 


In proceeding to the next load step, (N) , the temperature 
effects on the stress-strain curve are incorporated so as 
not to vialate this condition while holding the strains 
constant . 

'’9(n-l),(N) “If ° 


Thus, we are requiring that the change of temperature alone 
not effect the inelastic condition of the material. We 
accomplish this by realizing that 
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Therefore, by requiring that 


deP = 0 
Ij 


We force 

If-dT = 0 


(29) 
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This then gives us the mechanism for positioning our new yield 
surfaces in step, N. The step, N, solution then proceeds by applying 
the loads and boundary conditions and iterating to obtain 


dg 




within your specified convergence tolerance. 


CREEP ANALYSIS 


( 37 ) 


The creep analysis utilizes one of two possible creep representations. 
When tertiary creep is not considered to be of importance, the equation used 
is 


, 0 n^m ^ o r^ 
k e t + q e t 


( 38 ) 


where 

(T = a /lOOOOO, a = effective stress 
e e e 

k, m, n, q, r = material-dependent and temperature-dependent creep coef- 
f lei ents . 

Wien the material exhibits a significant amount of tertiary creep capa- 
bility, an alternate representation is used. Primary creep is represented by 
the Bailey-lloi ton law. 


P 

c 

c 


A2 A3 

Aj t 


( 39 ) 


Secondary creep is modeled with the expression proposed by Marin, Pao, and 
Cuff (Reference 19) 
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Tertiary creep is represented with an equation of the form 


T . ^9 - ^lOt 

“V “e 


(h\) 


^2» ••• Mo “ material-dependent and temperature dependent cri-> 
coefficients. 


CYANIDE also contains an orthotropic creep formulation. The creep s:*.*ain 
rate is assumed to be given by 


=ij “ ®ijkl®kl 


m 


where 

e-. “ strain rate tensor 
i-J 

“ stress tensor 


®ijkl 


Tensor whose components are functions of temperature, dg, and 
hardening rule and are derivable from input creep curves. 


The user can select from time hardening, strain hardening, or life frac- 
tion creep rule, depending upon the actual material characteristics. Strain 
hardening is ordinarily adequate for describing hardening behavior, providing 
that stress reversals do not occur. A stress reversal is considered to occur 
when 



(^♦ 3 ) 


Where eij® is creep strain measured from its current origin. When a reversal 
occurs, the origin is changed and the analysis proceeds (Reference 20} . 

The combination of general creep equations and creep rule makes the pro- 
gram very general In Its application to structures which undergo time-dependent 
plastic flow In which transient effects are not significant. 
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CYANIDE COMPUTER PROGRAM 


Many of the steps in the CYANIDE nonlinear finite element 
computer program are the same as those for a linear finite element 
analysis. The nonlinear effects are Introduced into the system of 
finite element equations by adding vectors of pseudoforces to the 
right hand side. 


|k| {6} = {F} + {F } + {F^} 

pc 

where 


|K| 

is 

{6} 

Is 

{F} 

Is 


Is 


Is 


the elastic stiffness matrix. 

the vector of nodal displacements. 

the force vector including thermal terms. 

the plastic pseudoforce vector. 

the creep pseudoforce vector. 


m 


For each Increment of loading, the nonlinear pseudoforces are iterated 
upon until the requirements of equl 1 ibri urn, compatabli ity, and the 
constitutive equations are met within user specified tolerances. Since 
this method does not require modification of the stiffness matrix during 
iterations It is very economical. This economy is magnified during 
cyclic analysis. The stiffness matrix need only be regenerated If the 
material properties are revised by thermal variation or if elements have 
been added or removed. 

MULT I AXIAL, VARIABLE TEMPERATURE EXAMPLE 

In a previous NASA contract, we Investigated one of the common thermal 
stress problems In AGTE's: turbine blade tip cracking. In that case, the 
critical region was shown by analysis and confirmatory testing tb have the 
cyclic stress-strain behavior noted in Figure 1. High temperature, time 
dependent flow rapidly relaxes the compressive stress such that on cool- 
down high tensile stresses are genrated. This process shakes down very 
rapidly to an almost elastic hysteresis loop based on modulus changes. 

In that case the problem was almost totally uniaxial in nature. 
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A second type of thermal stress problem prevalent in AGTE's is the hot 
spot. In this case, the stress strain response is definitely multiaxial. We 
will investigate a hot spot on a combustor shingle as being typical of these 
problems. Figure 2 shows a shingle segment. Taking advantage of its large 
radius of curvature and thinness, it was modeled as a flat plate in a 
condition of plane stress. The model is shown in Figure 3. Figures 4, 5 and 
6 show the nature of the hot spot at peak temperature and Figure 7 shows the 
heat-up cool-down temperature cycle at the center of the hot spot. This cycle 
was analyzed assuming no time dependent effects occcurred during heat-up and 
cool-down but that a one hour hold time was associated with the peak of the 
hot spot. 

The stress-strain results of the first cycle are shown in Figures 8, 9 
and 10 for the center of the hot spot. Figure 8 shows effective stress versus 
effective strain and Figures 9 and 10 show the biaxial stresses versus 
strains. Once again the effect of plasticity and creep is to generate tensile 
stresses during the cool-down portion of the cycle. The next series of 
figures shows the shakedown stress-strain results for the center of the hot 
spot. Figure 11 shows the effective stress versus effective strain shakedown 
values and Figures 12 and 13 show the shakedown biaxial stress cycle at the 
center of the hot spot. Thus this multiaxial thermal stress case, just as the 
uniaxial case, shakes down to almost elastic cycling with a high tensile mean 
stress. In addition, the stresses are almost proportional. These types of 
analyses are important in indicating the types of response and life tests 
needed . 
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• Arrows Denote Flxlt; In Indicated Direction 
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Figure 6. Temperature Contours At Peek Of Hot Spot 
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TEHPERATURE US. LOAD CASE (HSR802) CYANIDE ANALYSIS 
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Figure 9> Transverse Stress Versus Strain At The Center of the Hot 
Spot For the First Cycle 


Z-STRESS US. Z-STRAIN (HSR892) CVAHIDE ANALYSIS 
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Figure 10. Longitudinal Stress Versus Strain At the Center of The Hot 
Spot For the First Cycle 
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figure II. Shakedown Cycle of Effective Street Vertut Effective Strain At 
the Center of the Hot Spot 
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Figure 12. Shakedown Cycle For Trantvene Streit Venus Strain At the 
Canter of the Hot Spot 
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Z-STRCSS US. Z-STRAIN (HSR806) CVANICE ANALYSIS 



Figure 13. Shiktdown Cycle For Longitudinel Stress Versus Strein et the 
Center of the Hot Spot 


R-STRCSS US. 2-STRESS IHSRSSS) CYANIDE ANALYSIS 



Figure Ik. Shakedown Biaxial Stress Cycle At the Center Of the 
Hot Spot 
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